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ABSTRACT 

Atacama Pathfinder Experiment (APEX) mapping observations in CO (9-8) 
and (4-3) toward a high-mass star-forming region, NGC 6334 I, are presented. 
The CO (9-8) map has a 6."4 resolution, revealing a ~0.5 pc, jet-like, and bipolar 
outflow. This is the first map of a molecular outflow in a THz line. The CO (9-8) 
and (4-3) lines arising from the outflow lobes both show extremely high-velocity 
line wings, and their ratios indicate a gas temperature greater than 100 K and 
a density higher than 10 4 cm -3 . The spatial- velocity structure of the CO (9-8) 
data is typical of a bow-shock-driven flow, which is consistent with the associ- 
ation between the bipolar outflow and the infrared bow-shaped tips. In short, 
the observations unveil a highly-excited and collimated component in a bipolar 
outflow that is powered by a high-mass protostar, and provide insights into the 
driving mechanism of the outflow. Meanwhile, the observations demonstrate that 
high-quality mapping observations can be performed with the new THz receiver 
on APEX. 

Subject headings: ISM: jets and outflows — stars: formation — stars: early-type 



1. Introduction 

CO low-J transitions have been widely observed in studies of molecular clouds and 
bipolar outflows. These lines, e.g., CO (1-0) and (2-1), are easily excited at very low tem- 
peratures and mostly trace relatively cold molecular gas. Observations of high- J CO lines, 
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in particular those falling in the THz frequency regime (J up > 9), are very rare. These highly 
excited lines can be ob served with spaceborne observatories, e.g., In frared Space Observatory 
( Giannini et al. 20011) and most recently, Herschel Space Telescope (Ivan Kempen et al. 2010 : 
Yildiz et al. 20101 ) . On the other hand, ground-based THz observations, which can provide 
better angular resolutions down to a few arc seconds, are challenging even on extremely dry 
sites. To date only a very few ground telescopes have successfull y obtained astronomical 
THz CO spectra (e.g.. iKawamura et al. 20021 ; IWiedner et al. 20061 ) . However, all these ob- 
servations were carried out targeting one or a few points in an observed region, and had 
coarse spectral and/or spatial resolutions. Marrone et al. (2004) presented a CO (9-8) map 
of the OMC-1 region, but with a poor spatial resolution of 84". Therefore, morphological 
and kinematical information of the gas probed by THz CO lines is missing. 



As part of the science demonstration observations of the new THz receiver on APEXg, 
we mapped the CO (9-8) (1.037 THz) e mission in NGC 6334 I, the brightest far-IR source in a 
relatively nearby (1.7 kpc. lNeckel 19781 ) molecular cloud/Hn region compl ex (IMcBreen et al. 19791) 



NGC 6334 I has been extensively stud i ed over the last decade (e.g.. iHunter et al. 2 006 



Beuther et al. 20071 : iBeuther et al. 20081 : Emprechtineer et al. 2010l ; lvan der Wiel et al. 2010 . 
and references therein). It co ntains an ultracompac t (UC) HII region first detected in 6 



cm continuum (labeled "F" in iRodriguez et al. 19821 ) and appears as a "proto- Trapezium" 



system in 1.3 mm continuum (jHunter et al. 2o"o6 ). A bipolar molecular outflow, which is 



CO lines and in SiO and CS lines ( 


Bac 


filler & Cernicharo 1990: 


McCutcheon et al. 2000: 


Ridee & Moore 2001: 


Leurini et al. 2006 


). Here the CO (9-8) line was observed with a rel- 



atively high angular resolution. The line has an upper level energy (E up ) lying ~250 K 
above the ground and a critical density (n cr ) of order 10 6 cm -3 , thus is sensitive to hot and 
dense gas in shocked regions. In addition, the CO (4-3) line was simultaneously obtained 
with the dual-frequency receiver, enabling excitation analysis without significant effect from 
calibration uncertainties. 



2. Observations 

The observations were conducted on 2010 July 9, with ~0.25 mm precipitable water 
vapor, using a dual-frequency receiver system on APEX. The receiver system includes a 1 
THz channel, a copy of Herschel /HIFI band IV, and a co-aligned 460 GHz channel, designed 
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for pointing purpose (ILeinz et al. 2Q1CH ) . Two Fast Fourier Transform Spectrometers, based 
on the design presented by Klein et al. (2006), were configured to provide a spectral resolution 
of 122 kHz (~0.08 kms -1 ) over a 1.6 GHz bandwidth for the 460 GHz channel and a 
resolution of 183 kHz (~0.053 kms -1 ) over a 2.4 GHz bandwidth for the 1 THz channel. 
The system temperatures were around 360 K at 460 GHz and 2100 K at 1 THz. The 
beam efficiencies were determined to be 64% ± 8% and 21% ± 7% at 460 GHz and 1 THz, 
respectively, by observing Mars and Uranus. 

We performed on-the-fly (OTF) mapping observations, covering an area of about l'xO'5. 
Guided by the outflow elongation revealed by previous low- J CO observations, the long axis 
of the mapping area was tilted by 45° east of north. The effective integration time on 
each grid of the OTF map was 13 s. We pointed every hour toward (R.A., decl.)j 2 ooo = 
(1 7 h 20 m 53?44, -35°46 / 57. // 9), the peak position of the brightest 1.3 mm source (i.e., I-SMAl 



in 



Hunter et al. 20061 ) in the region, with the 460 GHz channel in the continuum mode. Fur- 



thermore, we constructed a pseudo 1 THz continuum map by averaging the line-free channels 
over an effective bandwidth of ~1.5 GHz, compared that pseudo continuum with a 345 GHz 
continuum map (convolved to the APEX beam of 6. "4) obtained from recent Submillimeter 
Array observations (Q Zhang, priv. comm.), and shifted the final maps by about 3."5 to 
match the high-accuracy 345 GHz peak at (R.A., decl.) J20 oo = (17 h 20 m 53?34, -35°46'59."1). 
Only CO (9-8) emission is clearly detected in the THz channel. In the 460 GHz channel, 
several lines are seen toward the 345 GHz continuum peak, but no line contamination is 
evident over the full observed velocity extent of the CO (4-3) emission. Unless specified, 
the data are presented in T£ scales. For a smoothed spectral resolution of 2 kms -1 , the rms 
noise levels are 0.27 K in CO (9-8) and 0.05 K in CO (4-3). 

Complementary IRAC 3-8 /iin images were retrieved from the Spitzer archive (Program 
ID: 30154). We use the Post Basic Calibrated Data products produced by the Spitzer Science 
Center pipeline S18.7.0. 



3. Results 

The CO (9-8) and (4-3) emission arising from the outflowing gas is detected with very 
high- velocity line wings, reaching LSR velocities —90 kms -1 in the blueshifted lobe and 60 
kms -1 in the redshifted lobe (the ambient cloud velocity is about —7.6 kms -1 ). Figured] 
shows the emission integrated in three velocity intervals, which are referred as the slow-wing 
(s-wing), fast- wing (f-wing), and extremely high- velocity (EHV) components following the 
terminology of Tafalla et al. (2010). Focusing on the high-resolution CO (9-8) maps, a 
bipolar outflow centered around the Trapezium-like millimeter continuum sources is clearly 
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detected. Provided a pointing uncertainty of 2-3", we refrain from an identification of the 
outflow driving source from the millimeter sources. The outflow extends about 0.5 pc from 
the northeast (NE) to southwest (SW), in an orientation consistent with low- to mid- J CO 
observations. A blueshifted structure extending from the millimeter sources to the northwest 
also shows EHV emission, and is likely attributed to another outflow. Since the feature 
extends beyond our mapping area, it is not further analyzed in this work. 

Figures [2ja), [2(b) show the CO (9-8) and (4-3) spectra extracted toward the outer 
peak on each of the blue- and redshifted lobes (marked as "B" and "R" in Figure [T]) . The 
two transitions have similar profiles: the s-wing component has its brightness fast decreasing 
as the velocity increases; the f-wing component appears as a plateau with relatively flat 
brightness variation; the EHV component has greater CO (9-8)/(4-3) ratios, and shows a 
secondary peak in the blueshifted lobe. The CO (9-8)/(4-3) ratios averaged within each 
of the three velocity intervals are measured to be 0.70 (s-wing), 0.69 (f-wing), 0.92 (EHV) 
toward "B" and 0.79 (s-wing), 0.81 (f-wing), 0.92 (EHV) toward "R". Again, greater ratios 
are found in the EHV component, which imply a higher excitation condition for the EHV gas. 
We p erformed large- velocity-gr adient (LVG) calculations on the line ratios using the RADEX 



code (Ivan der Tak et al. 20071 ). Figure [2(c) shows the CO (9-8)/ (4-3) ratios as a function 
of gas density and kinetic temperature. The adopted CO column density to line width ratio 
is 10 16 cm _2 /(kms _1 ), corresponding to a total column density on the order of 10 17 cm -2 as 
constrained by Leurini et al. (2006) based on 13 CO observations. The calculations show that 
for the observed range of line ratios, the gas density is > 10 4 cm -3 and the temperature is 
> 100 K. To further constrain the physical conditions of the outflowing gas, it is instructive 
to adopt a typical CO abundance of 10~ 4 and assume that the gas dimension along the 
light-of-sight is similar to that in the plane-of-sky. In this situation the CO column density 
and H 2 density are connected by a size scale, and that scale is used to estimate the beam 
filling factor. We then ran a series of LVG models exploring the parameter space of kinetic 
temperature, H 2 density, and CO column density, and compared the brightness temperatures 
predicted by the models with the corrected antenna temperatures (T m b divided by a beam 
filling factor). We find that the models with a density of ~2-3 x 10 4 cm -3 and a kinetic 
temperature of ~400 K appear to reasonably agree with the observed EHV emission, while 
the models with a similar density but a lower temperature of ~300 K are able to reproduce 
the s-wing/f-wing emission. 

Figure [3] shows a comparison between the CO (9-8) emission and shock activities 
probed by mid-IR and NH 3 maser observations. In Figure E](a), the Spitzer three-color 
composite image reveals a pair of bow-shaped tips around the CO (9-8) peaks in ex- 
cess 4.5 u m (green) emission, which often traces jets and bow-shocks in protostellar out- 



flows (e.g.. iNoriega-Crespo et al. 20041 ; ISmith et al. 20061 ; lOiu et al. 20081 ). Several authors 
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have detected 2.12 / im H2 knots in this region ( IDavis fc Eisloffel 19951 ; iPersi et al. 1996 



Eisloffel et al. 20001 ; ISeifahrt et al. 20081). Two bow-sh aped H 2 knots along the CO out- 



flow axis (denoted "a" and "/?" in ISeifahrt et al. 20081 ) coincide with the 4.5 /im emission 
knots reported here . This suggests that the excess 4.5 nm emission can be mostly due 
to shocked H 2 lines (ISmith fc Rosen 20051 ; ISmith et al. 20061 ; iDe Buizer fc Vacca 20101 ). but 
the contribution from CO v=l-0 or Hi Bra lines cannot be ruled out without spectroscopic 
observations. Furthermore, Figure Qfb) shows the 4.5 /im to 3.6 urn , emission ratio to high- 
light the shocked emission (ITeixeira et al. 20081 ; iTakami et al. 20101 ). While field stars have 
the brighter emission falling in the 3.6 /im band, shocked emission appears most prominent 
in the 4.5 /im band, thus Figure 0(b) not only shows the bow-shaped tips seen in Figure 
[3ta) but also unveils a new feature which is suppressed in Figure 0(a) by a nearby bright 
star. Each of the th ree 4.5 /im excess features is associated with an NH 3 (3,3) maser spot 
(IBeuther et al. 20071 ). Compared to the 1st moment (intensity- weighted velocity) map of 
the CO (9-8) emission, the shocked features seen in IR and maser emission coincide with 
the outer parts of the outflow lobes where most of the gas moves at high velocities. 



4. Discussion 

The most surprising result of our observations is the detection of high- velocity CO (9- 
8) emission (f-wing and EHV) over a spatial extent of order 0.5 pc. This provides direct 
evidence for the presence of a highly-excited component in a bipolar outflow. The dynamical 
timescale of the CO (9-8) outflow, deriv ed from the flow radius of ~0.22 pc divided by the 



outflow terminal velocity of ~77 kms 1 (ILada 19851 ). is about 2.6 x 10 3 yr. This agrees with 



previous estimates of 2-4 x 10 3 yr based on the CO (2-1), (3-2), and (4 -3) observations 



(IBachiller fc Cernicharo 1990c iMcCutcheon et al. 2000l ; iLeurini et al. 20061 ). Therefore the 



outflow is very young, preferring an evolutionary stage prior to the formation of an UC 
Hn region for its driving source. Compared to the mid-IR and NH3 maser observations, 
the outflow seems to be closely related to bow shocks. The CO (9-8) observations have a 
~ 6" resolution, allowing us to further explore the kinematics of the outflow and its driving 
mechanism. 

Figure 0] shows the position-velocity (P-V) diagram of the CO (9-8) emission, con- 
structed approximately along the outflow axis. Overall, the outflow velocities (relative to 
the ambient cloud velocity) are greater at positions farther away from the central source. 
Such a velocity structure is also discernable in the lst-moment map shown in Figure O^b). 
In particular, the maximum velocities in the redshifted lobe roughly linearly increase with 
the distances from the central source, following a "Hubble" law. Theoretically, outflows 
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forming from jet-driven bow shocks have been shown to exhibit Hubble-like P-V relations, 
where the f orward velocities of the gas around the bow's head are great e r than along th e 
bow wings Jzhang fc Zheng 19971 : ISmith et al. 19971 bownes fc Rav 19991 : iLee et al. 200lh . 
On the other hand, Stahler (1994) obtained the Hubble P-V relation in a model of ambi- 
ent gas entrained by a turbulent jet. However, as pointed out by Lada & Fich (1996), a 
critical assumption of Stahler's calculation is that the gas mass in every slab perpendicular 
to the jet axis has the same power-law distribution with velocity, which is not seen in our 
CO (9-8) or (4-3) data. The P-V pattern we observe is very similar to that reproduced 
by the bow-shock calculations of Zhang & Zheng (1997) and simulations of Smith et al. 
(1997) and Downes & Ray (1999). We also note that the P-V pattern shown in Figure 
H] closely resembl es that of the NGC 2264G outflow, wh ich is a Class outflow driven by 
a precessing jet (ILada &: Fich 19961 : iTeixeira et al. 20081 ). Moreover, the bow-shaped tips 
seen in the IR independently suggest the association between bow shocks and the CO (9-8) 
emission. Thus, we suggest that the gas probed by the CO (9-8) line wings is accelerated 
by bow shocks created by a jet impacting the ambient material. The nature of the EHV 
gas, however, is less clear. From Figure |2j the terminal velocities of the EHV component 
reach > 100 kms -1 if the projection effect is taken into account. Non-dissociative C-type 
shocks are expected to accelerate ambient gas to velocities up to ~ 50 kms -1 . A recent 
chemical investigation of EHV low-mass outflows chal lenges the hypothes is that molecules 
in EHV gas reform behind dissociative J-type shocks (ITafalla et al. 20101 ). It is likely that 
the EHV gas has an origin close to the centr al source rather than being accelerated ambient 
gas jQiu fc Zhang 2009t KafaJla et al. 201~oL 



Most low- J CO observations of molecular outflows show power-law intensity- velocity 
(I-V) relations, Ico(v)ocv ~ 7 , with the power-law index, 7, typically around 2 and in some 



cases steepening to larger values at velocities higher than 1 0-30 kms 1 (e.g., IStahler 1994 



Lada fc Fich 19961 : iDavis et al. 19981 : iRidge fc Moore 20011 ). If the line wings are optically 



thin and the CO abundance and excitation temperature do not vary much, as assumed 
in many studie s, the shape of I-V plots represents mass- velocity distributions (see, e.g., a 
compilation by lRicher et al. 20001 ) . Th e I-V property provides an important test for mo dels 



of molecular outflow acceleration (e.g.. iDownes fc Cabrit 2003c iKeegan fc Downes 20051 ). In 



Figure EJ we plot the spatially integrated CO (9-8) intensity as a function of velocity. For 
v < 55 kms -1 , the intensity distribution can be described by power-laws with 7 = 0.78±0.04 
and 0.90 ± 0.04 for the blue- and redshifted lobes, respectively. At higher velocities, the 
intensity rapidly drops, and the power-law fits give 7 = 2.59 ± 0.33 (blue) and 4.32 ± 1.27 
(red). The I-V relation in CO (4-3) is similar, with 7 = 1 in each lobe up to ~ 55 kms -1 
and steepening further out. In addition, Ridge & Moore (2001) reported 7 « 1 for v < 50 
kms -1 with the CO (2-1) observations of this outflow. Thus, the I-V relation of the NGC 
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6334 I outflow, measured in low- to high- J CO lines for velocities up to ~ 55 kms -1 , appears 
to be among the shallowest known and is shallower than most theoretical predictions (see a 
recent review by lArce et al. 20071 ). Simulations of collimated o utflows driven by den se jets 
predict 7 close to 1 (jet to ambient gas density ratio set to 10, iRosen Smith 20041 ). But 
the thrust of the simulated jets can only drive low-mass outflows, and the 7 ~ 1 power-law 
in the simulations occurs at velocities up to 10 kms -1 , far lower than the value observed 
here (55 kms -1 ). It is therefore yet to be confirmed whether the observed shallowness in the 
I-V relation can be a consequence of the high density of an underlying jet. Finally, Ico( v ) a t 
> 55 kms -1 behaves dramatically different from that at lower velocities. It probably reflects 
a different nature of the EHV g clS clS discussed above. 



5. Summary 

The observations reveal a jet-like and extremely high-velocity outflow originating from a 
high-mass protostar. The CO (9-8) map presented here is the very first THz frequency map 
of a molecular outflow. It indicates that a highly-excited CO line (E np ~ 250 K, n CT ~ 10 6 
cm -3 ) can still have broad line wings tracing parsec scale outflows. LVG calculations of the 
CO (9-8)/(4-3) line ratios give kinetic temperatures > 100 K and H 2 densities > 10 4 cm -3 . 
The spatial-velocity structure of this relatively hot and dense outflowing gas suggests that 
the outflow is driven by jet bow-shocks. The detection of IR bow-shaped tips coincident 
with the CO outflow lobes further supports this scenario. It is worth noting that the outflow 
shows an exceptionally shallow I-V relation in low- to high- J CO lines, though whether 
that shallowness is caused by a high density of the jet is to be addressed. The outflow also 
possesses EHV gas. The CO (9-8) / (4-3) ratios and the I-V relation of the EHV emission 
both appear to be different from lower velocity line wings, which may suggest a different 
origin (close to the central protostar) of the EHV gas. 

Finally, the observations demonstrate that under excellent weather conditions, high- 
quality THz observations can be performed with the new THz receiver on APEX. 



This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with the National Aeronautics and Space Administration (NASA). 
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Fig. 1. — APEX CO (9-8) emission integrated in three velocity intervals: (a) s-wing, the 
lowest contour is 15.6/18.0 Kkms -1 (blue/red) and the contour spacing is 3.9/4.5 Kkms -1 
(blue/red); (b) f-wing, the lowest contour is 16.5/14.5 Kkms -1 (blue/red) and the contour 
spacing is 6.6/5.8 Kkms -1 (blue/red); (c) EHV, the lowest contour is 11.6/9.6 Kkms -1 
(blue/red) and the contour spacing is 5.8/4.8 Kkms -1 (blue/red), (d)-(f) The CO (4- 
3) emission integrated in the same velocity intervals; the equally spaced contour levels go 
from 25% to 95% of the peak intensities, which are 173/158 Kkms -1 (blue/red) for s-wing, 
182/137 Kkms -1 (blue/red) for f-wing, and 76/38 Kkms -1 (blue/red) for EHV. Four filled 
squares denote millimeter continuum sources detected by Hunter et al. (2006). A small 
circle labels the shifted THz continuum peak. Two "+" symbols mark approximately the 
outer peaks of the blue- and redshifted lobes, namely "B" and "R", respectively. A solid 
line in panel (a) shows the cut of the position- velocity diagram (see Figure @J. 
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Fig. 2. — (a) The CO (9-8) and (4-3) spectra extracted from the "B" position (Figure [Q 
shown in gray and black histograms, respectively; the CO (9-8) data have been convolved to 
the CO (4-3) beam, which is 14."4 at FWHM; the vertical dashed lines delineate the three 
velocity intervals (see Figure [T] and Section [3]). (b) Same as (a), but extracted from the 
"R" position (Figured]), (c) Curves of constant CO (9-8) /(4-3) ratios as a function of H 2 
density and kinetic temperature, for an adopted CO column density to line width ratio of 
10 16 cm- 2 /(kms- 1 ). 
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Right Ascension (J2000) 

Fig. 3. — (a) Spitzer three-color composite image with the 3.6, 4.5, and 5.8 /im emission 
coded in blue, green, and red, respectively; dotted lines outline two areas where the 4.5 
/xm/3.6 /im band ratio is computed; other symbols are the same as shown in Figure [TJ 
(b) The colored image shows the 1st moment map of the CO (9-8) emission, and the gray 
contours show the Oth moment map with the lowest contour of 110 Kkms -1 and the contour 
spacing of 22 Kkms -1 ; the gray scale represent the 4.5 /iin/3.6 /im band ratio at levels of 
1.7, 1.85, 2.0, and 2.15; white contours represent the NH 3 (3,3) maser emission originally 
presented in Beuther et al. (2007). 
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Fig. 4. — CO (9-8) position-velocity diagram constructed along the direction going through 
the peaks of the outflow lobes (see Figure (Ha)); the lowest contour and the contour spacing 
are 0.75 K; the offsets are measured with respective to the midpoint of the two peaks, and 
should roughly represent the distances to the outflow central source. 
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Fig. 5. — Spatially integrated intensity of the CO (9-8) emission as a function of outflow 
velocity. Filled circles and open squares denote the blue- and redshifted emission measured 
within each 2 km s" 1 velocity channel, respectively. Power-law fittings to the data are overlaid 
in solid (blue) and dashed (red) lines. 



